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Abstract 

This paper deals with the mathematical modelling of large strain magneto¬ 
viscoelastic deformations. Energy dissipation is assumed to occur both due to the 
mechanical viscoelastic effects as well as the resistance offered by the material to 
magnetisation. Existence of internal damping mechanisms in the body is consid¬ 
ered by decomposing the deformation gradient and the magnetic induction into 
‘elastic’ and ‘viscous’ parts. Constitutive laws for material behaviour and evolution 
equations for the non-equilibrium fields are derived that agree with the laws of ther¬ 
modynamics. To illustrate the theory the problems of stress relaxation, magnetic 
field relaxation, time dependent magnetic induction and strain are formulated and 
solved for a specific form of the constitutive law. The results, that show the effect 
of several modelling parameters on the deformation and magnetisation process, are 
illustrated graphically. 
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1 Introduction 


Magnetorheological elastomers (MREs) are materials that change their mechanical be- 
havionr in response to the application of an external magnetic held. These elastomers 
have received considerable attention in recent years dne to their potential nses as variable 
stiffness actnators for mechanical systems with electronic controls. MREs are particn- 
larly nsefnl for their tnneable elastic modnlus and a rapid response to the magnetic held, 


cf. Bose et ah (2012). A common preparation method is mixing magnetically permeable 


particles into liqnid monomer and letting the mixtnre to polymerise. Curing, when done 
in the absence of magnetic held, results in an isotropic material while curing in the pres¬ 
ence of a magnetic held causes the particles to align in a particular direction and results 
in a material with a directional anisotropy. The ferromagnetic particles are usually be¬ 
tween 1-5 /im in size and kept between 0-30% by volume of the entire mixture. Such 


elastomers have been reported to be prepared and analysed by Jolly et ah (1996), Cinder 
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et al. 

(2002 , 

Varga et al. 

(2006) 

Boczkowska and Awietjan 

(^2009 

), and 

Bose and Roder 


( |2009D . 

Mathematical modelling of the coupling of electromagnetic helds in deformable ma¬ 
terials has beeen an area of active research in the past, see, for example, the works of 
Pao (1978) and Eringen and Maugin (1990). Recently, a new constitutive formulation 


based on a ‘total’ energy density function has been developed by [Dorfmann and Ogden] 
(2003, 2004), wherein the solutions of some boundary value problems were obtained us¬ 
ing different energy densities. It has been shown that any one of the magnetic induction 
vector, magnetic field vector, or the magnetisation vector can be used as an independent 
variable of the problem and the other two obtained through the constitutive relations. 
The relevant equations used by them are based on the classic work of Pao (1978) in which 
the equations of motion for an isotropic non-polar continuum in an electromagnetic field 
are described by Maxwell’s equations and the mechanical and thermodynamical balance 
laws. This formulation has been particularly useful in recent years in dealing with prob¬ 
lems related to magnetoelasticity and using this, further boundary value problems on 
nonlinear deformation and wave propagation have been studied by [Bustamante et ah 
(2007), Ottenio et al. (2008) and Saxena and Ogden (2011). 

The response to an applied magnetic induction is, however, not exactly instanta¬ 
neous for all materials. On the application of a sudden external magnetic induction, the 
magnetic field (or equivalently the magnetisation) developed inside the material is not 
constant. Starting with some initial non-equilibrium value, it gradually approaches equi¬ 
librium in some finite time (say ti) depending on the existing deformation and various 
material parameters. The synthetically developed magnetoelastic materials are usually 
polymer based, hence also viscoelastic in nature. Thus there is development of a viscous 
overstress on deformation or on the application of a body force that vanishes after a time 
^2 which is usually different from ti above. This time-delay in response is a very important 
factor to consider while designing electromechanical actuators from magnetorheological 
elastomers. Thus, these two forms of dissiption - due to mechanical and due to magnetic 
effects need to be modelled appropriately. In order to consider the magnetic and mechan¬ 
ical dissipation effects, the previously stated theory of magnetoelasticity by [Dorfmann 


and Ogden (2004) is generalised by combining with the existing theory of mechanical 


viscoelasticity. 

Viscoelastic material modelling can generally be classified into two main classes, i.e. 
purely phenomenologically-motivated and micromechanical based network models. In 
some literature, the viscoelastic modelling approach is also divided, on the one hand, due 
to the nature of the time-dependent part of the stress, on the other hand, due to the 
nature of the evolution equation. The phenomenological modelling approach can also be 
distinguished based on the type of internal variables, i.e. stress type internal variables in 
the form of convolution integrals, cf. Simo (1987), Holzapfel and Simo (1996), Lion (1997), 


Kaliske and Rothert (1997); and strain-type internal variables that originate from a 


multiplicative decomposition of the deformation gradient, cf. Reese and Govindjee ( 1998[ ) 
and Huber and Tsakmakis (2000). The latter group decomposes the deformation gradient 
into elastic and inelastic parts where the inelastic part is determined from a differential 
type flow rule. In both the cases, the total stress is decomposed into a viscosity induced 
overstress and an equilibrium stress that corresponds to stress response at an inhnitely 
slow rate of deformation or the stress response when the time-dependent viscous effects 
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are completely diminished. Within the setting of the multiplicative decomposition of the 
deformation gradient, Reese and Govindjee (1998) proposed an evolution law which, when 
linearized around the thermodynamical equilibrium, yields the hnite linear viscoelasatic 
model of Lubliner] (1985). Koprowski-Theiss et ah (2011) proposed a nonlinear evolution 
law, which after being proved to be thermodynamically consistent, has been used in this 
paper. For the modelling based on stress-type internal variables, the time-dependent 
overstress part is expressed as an integral over the deformation history, cf., Simo (1987) 
and Amin et ah] (2006). 

The second class of viscoelastic material modelling is based on micromechanical the¬ 
ories derived using the underlying molecular structures, see, for example, the works of 


Bird et al. 

(] 

.987 

), 

Doi and Edwards 

(1988 

), 

Bergstrom and Boyce 

(1998 

) and 

Miehe 

and Goktepe 

(2005 

1 . These have been developed over the years to describe the viscous 


behaviour of molten polymers and physically cross-linked rubber-like materials. The 
bead-spring model of Bird et al. (1987), the reptation-type tube models of |de Gennes 


(1971) and Doi and Edwards (1988), and the transient network models of Green and 


Tobolsky (1946) can be mentioned as examples in this area. The theory for transient 


network models explains the stress relaxation phenomena as a consequence of breakage 
and reformulation of the polymer cross-links constantly, cf. Green and Tobolsky ( 1946| 
and Reese (2003). Reptation-type tube models are developed for the dehnition of the 
motion of a single chain in a polymer gel. The constraints on the free motion of a single 
chain are qualitatively modelled as a tube-like constraint and the motion of the chain is 
described as a combination of Brownian motion within and reptational motion along the 
tube. Recently, a growing interest can be observed to combine these approaches which 


yield the so-called micromechanically motivated models, see, for example, Linder et al. 

ponj ). 

As a hrst step in the magneto-viscoelastic modelling, we model an isotropic material 
and take a phenomenological approach based on the multiplicative decomposition of the 
deformation gradient in line with Lubliner (1985). An additive decomposition of the mag¬ 
netic induction vector into equilibrium and non-equilibrium parts is proposed to model 
magnetic dissipation phenomena. The equilibrium part of the energy is taken to be a 
generalisation of the Mooney-Rivlin elastic model to include magnetic effects, while the 
non-equilibrium part is a slightly simplihed version that looks like a neo-Hookean type 
magnetoelastic model. Using a Glausius-Duhem form of the second law of thermody¬ 
namics, we obtain evolution equations for these physical quantities to be able to perform 
numerical calculations. 

This paper is organised as follows. In Section 2, the theory of nonlinear magneto¬ 
viscoelasticity is presented taking into account the case of hnite deformation. Starting 
with the governing Maxwell’s equations and the laws of momentum balance, we show the 
existence of a total stress tensor of Dorfmann and Ogden (2003). The deformation gra¬ 
dient and the magnetic induction are decomposed into equilibrium and non-equilibrium 
parts. Using the laws of thermodynamics and a form of the Helmholtz free energy func¬ 
tion, constitutive equations are derived along with the conditions to be satished by the 
evolution equations of the non-equilibrium quantities. 

In Section 3, for the purpose of obtaining numerical solutions the energy density 
function and the evolution equations for the non-equilibrium quantities are specialised 
to specihc forms. Several material parameters to model magneto-viscoelastic coupling 
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Fig. 1. The material, intermediate, and spatial configurations with the corresponding 
magnetic vectors and deformation tensors. 


are introduced in this step. In Section 4, we consider four different types of deformation 
and magnetisation processes to study the effects of the underlying magnetic induction, 
deformation, strain rate and magnetic induction rate on the total Cauchy stress and 
the magnetic field relaxation process. It is observed that changing the newly defined 
magneto-viscoelastic coupling parameters can affect the magnitude of the overstress, the 
excess magnetic field and their decay times. Initial deformation can affect the decay 
time and magnitude of the induced magnetic field. Effects of the deformation, applied 
magnetic induction and the material parameters on the computed physical quantities 
(such as stress and magnetic field) are illustrated graphically. Section 5 contains some 
brief concluding remarks. 


2 Theory of nonlinear magneto-viscoelasticity 


We consider an incompressible magnetoelastic material which, when undeformed and un¬ 
stressed and in the absence of magnetic fields, occupies the material configuration Bo with 
boundary dBo. It is then subjected to a static deformation due to the combined action 
of a magnetic field and mechanical surface and body forces. The spatial configuration at 
time t is denoted by Bt with a boundary dBt- The two configurations are related by a 
deformation function x which maps every point X G So to a point x = x(X, t) G Bt- 
The deformation gradient is defined as F = Grad x, where Grad is the gradient operator 
with respect to X. Its determinant is given by J = det F = 1 for the present case of 
incompressibility. 

To take into account mechanical viscous effects, we assume the existence of an in¬ 
termediate configuration Bi that is related to Bt by a purely elastic deformation and is 
related to Bq by a pure viscous motion. The intermediate configuration Bi is postulated 
only to model the dissipation effects. This is in parallel to the energy-conserving magne¬ 
toelastic deformation from Bq to Bt- Following Lubliner (1985) and Reese and Govindjee 


(1998), this motivates the decomposition of the deformation gradient into an elastic and 

F = F.F,. (1) 


a viscous part as 
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For future use we define the right and the left Cauchy-Green strain tensors as C = F*F 
and b = FF*, respectively. Similar quantities corresponding to F.^ and Fe are dehned in 
Bq, Bi and Bt as shown in Fig. [T] 

It is further assumed that the material is electrically non-conducting and there are no 
electric helds. Let cr be the ‘mechanical’ Cauchy stress tensor and r the total Cauchy 


stress tensor (see, for example, Dorfmann and Ogden (2004) for its dehnition), p the mass 
density, the mechanical body force per unit mass, a the acceleration of a point, ff the 
electromagnetic body force per unit volume, Ih the magnetic held vector, b the magnetic 
induction vector, and m the magnetisation vector. The following balance laws need to 
be satished 


div cr -|- ff -f p^rn = pa, div r -|- pf^ = pa. 


T = T, 


( 2 ) 


curl]h = 0, divb = 0. (3) 

In Eq. (|^, the hrst two equations are equivalent forms of the balance of linear momentum 
and the third is the angular momentum balance equation. Eq. (3)i is the specialisation 
of Ampere’s law appropriate to the present situation and Eq. (3)2 is the statement of 
impossibility of the existence of magnetic monopoles. Here and henceforth, grad, div, 
curl denote the standard differential operators in Bt while Grad, Div, Curl denote the 
corresponding operators in Bq. The three magnetic vectors are connected through the 
standard relation 

b = po [li m], (4) 


Po being the magnetic permeability of vacuum. The connection between cr and r is 


T = cr 




bob 


- [b ■ b] i 


-|- [m ■ b] i — b 


m, 


(5) 


where i is the second order identity tensor in Bt and we h ave used th e expression for the 
magnetic body force as ff = [gradb]*m; see, for example, Pao (1978). 

The total Piola-Kirchhoff stress and the Lagrangian forms of ffi, m, and b for an 
incompressible material (J = 1) are dehned by 


s = f-Vf-* 

= f-Vf-' + po^ 


- - [[CB] ■ B] C-^ 


+ [M ■ B] - 


I]C-S (6) 


H = F^li, M = F*m, B = F^^b. (7) 

We use the above relations to rewrite the governing equations ([^ and (|^ in terms of the 
Lagrangian variables as 


Div (SF*) -|- pfm = pa, S* = S, CurllH = 0, DivB = 0, 
while the relation (|^ becomes 

CB = Po [e + M]. 


( 8 ) 


(9) 


In magnetorheological elastomers, in addition to the mechanical viscoelastic effects, we 
propose that energy dissipation also occurs due to the resistance offered to magnetisation 
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of the material. On the sudden application of a constant magnetic induction, the magnetic 
held generated inside the material starts from an initial non-equilibrium value and then 
evolves to approach an equilibrium value. To model these effects, we assume the existence 
of a dissipation mechanism by including magnetic induction like ‘elastic’ and ‘viscous’ 
internal variables be and b„, respectively. Their Lagrangian counterparts are given by 
Be and such that 

b = be + b^, B = Be-|-B^. (10) 

The above additive decomposition of the magnetic induction is motivated by a similar 
decoupling into elastic and viscous parts of the deformation in viscoelasticity theory. Since 
magnetic induction is a vector, an additive decomposition is applied here as opposed to 
the multiplicative decomposition of the deformation gradient in Eq. ([^. The behaviour 
of the internal variables dehned in the above is assumed such that if a constant magnetic 
induction B is applied at time t = 0, then at that instant Be = B and B„ = 0 . As time 
progress, the magnetic induction is gradually and entirely transferred to B^. Thus, 

Be = B, B„ = 0 at f = 0, 

Bg — y 0 , Bj; —y B as t —y oo. (11) 


2.1 Thermodynamics and constitutive modelling 


In this section, starting with the laws of thermodynamics, we use the balance equations 
(|^ and (|^ to obtain constitutive laws that dehne the material behaviour. Necessary 
conditions of energy dissipation associated with the viscous and the magnetic dissipation 
processes are obtained that need to be satisfied by any magnetoelastic deformation to be 
thermodynamically admissible. 


Balance of energy is written in the local form as (cf. Pao (1978), Dorfmann and Ogden 


(2003)) 


(u + ) -h divQ = div (crv) -h [pf^ + f] - v + pR + Wg 


( 12 ) 


where U and R denote the internal energy and radiant heating per unit mass, Q is the 
heat flux, cr is the purely mechanical Cauchy stress, is the mechanical body force 
(assumed to be zero in the analysis later), and We is the electromagnetic power given for 
the present case as Wg = —m ■ dh/dt. 

Let S be the specihc entropy and "d be the temperature. On introducing a specihc 
Helmholtz free energy T through 


^ = U-^S, 


and using the Clausius-Duhem form of the second law of Thermodynamics 


dS , 




> 0 , 


we arrive at the following inequality 


d^ 


dF dh 

— -m ■ — 

dt dt 


> 0 . 


( 13 ) 


( 14 ) 


( 15 ) 


6 








The symbol : denotes a double contraction operation between two second order tensors. 
In the calculations above use has been made of Eq. and the temperature is assumed 
to be constant. 


We now introduce a total energy function similar to the one used by Dorfmann and 


Ogden (2004) 


r2(F, Cy, B, B^,) = pT(C, C.,;, FB, FB.^) + -—B • [CB]. 

2po 


(16) 


This considers the magnetic induction vector B as an independent quantity and leaves 
the magnetic held El to be determined using a constitutive law. The magnetisation M, if 
required, can be obtained using the relation ([^. 

The above equation, on differentiation with respect to time gives 




dQ 


^ dt dt 


—] 
2po 


while using ([^, we obtain 


dh 


— m ■ ^ = —i 


dt 


dC^ 

dt 


F-i—I 
dt 


- [«i ■ 

po dt 


dt 


Note that the constraint of incompressibility can be expressed as 


dt 


= C 


-1 


dt 


(17) 


(18) 


(19) 


For an incompressible material, dC/dt is not arbitrary, but the inequality (15) must 
be satished for every dC/dt governed by the above constraint. Consequently adding a 
scalar multiple of this zero term and substituting Eqs. ( [TT] ) and (18) to the inequality 
(15), a form of the second law of thermodynamics is obtained in terms of the ‘total’ Piola- 


Kirchhoff stress tensor and the physical quantities dehned in their Lagrangian description 
as 

- — + - S + ^ + H ■ — > 0, 

dt 2 ^ ^ dt dt 


( 20 ) 


p being a Lagrange multiplier associated with the constraint. It is noted here that one 
can equivalently use a nominal stress tensor T = F^^r instead of S and take F 


as 


an independent variable instead of C and substitute in (15). This would result in an 


inequality that yields the constitutive relations used by Dorfmann and Ogden (2004). 


Taking partial derivatives of D with respect to its arguments and substituting in the 
inequality above, we get 


on 


dC 

dt 

dC, 


+ 


~ 

90 dM,. 


9C, 


dt 9B,, dt 


dt 

> 0 , 


( 21 ) 


From the arguments of Coleman and Noll (1963), the following constitutive equations 
are obtained 

S = 2gf-pC-‘, H=|l, (22) 


dC 


9B’ 
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(23) 


along with the dissipation condition 

dQ dCy dfl dEy 

-: - - H- - < 0. 

dCy dt dt 


For the sake of simplicity in the computations later, we further assume that the 
non-equilibrium magnetic induction and the non-equilibrium strain tensor Cy are 
independent from each other. This reduces the above inequality to the following separate 
conditions 


dn dCy dn dMy 

-: < 0 ,- —- < 0 , 

dCy dt 9B^ dt 


(24) 


which should be satished by any magnetoelastic deformation process to be thermody¬ 
namically admissible. 

The total energy stored in the body can be split into an equilibrium part associated 
with the direct deformation from Bq to Bt, and a viscous part due to the internal variable 
Bg and the elastic deformation from Bi to Bt- This is a slightly general form of the purely 
mechanical energy decomposition by [Reese and Govindje 


lee 


(1998) 


n(C,C„,B,B^) 


r2g(C, B) -|- (C, C^,, B, B^). 


(25) 


Here, the viscous part of the energy depends on the viscous parts of the deformation and 
the magnetic induction. Thus the arguments of Qy can be equivalently changed as either 
one of ily (C, C^,Bg), ily (Cg,B,B^) or ily (Cg,Bg). 

Substituting this form of hi into the inequalities ( |24j ) we obtain the dissipation condi¬ 
tions that should be necessarily met in order to satisfy the second law of thermodynamics 


dQy _ dCy ^ 

dCy ''dt - ' 

d^ly (iB„ 

—^-^ < 0 . 

5B„ dt 


(26) 

(27) 


It is noted here that the above theory can easily be generalised to include multiple 
dissipation mechanisms in the body. In the case of M mechanical and N magnetic 


mechanisms, we may dehne F^, 

■pM. pi 

pM. 

V ) 


...,Bf; Bi,., 

such that 


■p _ 

^ ^ v") 

Vi = 

1,.., 

.,M, 

(28) 

B 


Vj = 

:1,.. 

.,N. 

(29) 


The dissipation condition to be satished in this general case is 


E 


dn 


dC\ 

dt 


N 

E 

i=i 


dn 

dni 


dEj 

dt 


< 0 . 


(30) 


3 Specialised constitutive laws 


With a motivation of obtaining numerical solutions to some magneto-viscoelastic defor¬ 


mation problems, we specialise the energy in (25) to specihc forms in this section. Evolu¬ 
tion equations for and B^ are also derived that satisfy the thermodynamic constraints 


(26) and (27). 
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3.1 Energy functions 

The material is assumed to be isotropic following which the equilibrium part of the 
energy density function is considered to be a generalisation of the classical Mooney- 
Rivlin function to magnetoelasticity of the form 


fie 


he 

T 


1 + Oe tanh 



[[1 + n] [R-3] + [l-n] [/ 2 - 3 ]] 

+qh + rid, 


(31) 


where /i,/ 2,/4 and Iq are the standard scalar invariants in magnetoelasticity (see, for 
example, Dorfmann and Ogden (2004|)) dehned as 


/i = C : I, 




/4 = |B ® B] : I, 
h = [ [CB] 0 [CB] ] : I, 


(32) 


I being the second order identity tensor in Bq- 

This is a slight generalisation of a Mooney-Rivlin type magnetoelastic energy function 
proposed by Ottenio et ah (2008). Here /ig is the shear modulus of the material in the 


absence of a magnetic field and n is a dimensionless parameter restricted to the range 
— 1 < n < 1, as for the classical Mooney-Rivlin model. The term [1 + CKetanh (/ 4 /me)] 
corresponds to an increase in the stiffness due to magnetisation and the phenomenon of 
magnetic saturation after a critical value of magnetisation. The parameter rUe is required 
for the purpose of non-dimensionalisation while Oe is a dimensionless positive parameter 
for scaling. The magnetoelastic coupling parameters q and r have the dimensions of /Tq 
F or = q = r = 0, this simplifies to the classical Mooney-Rivlin elastic energy density 
function widely used to model elastomers. 

Let the natural basis vectors in Bq be identihed with a set of covariant basis vectors 
{Ga} and its dual basis with a set of contravariant basis vectors {G“}, a e {1,2,3}. 
Similarly defining {ga} and (g"} for Bt and {ga} and (g"} for Bi, we obtain the following 
forms for the deformation gradient and the right Cauchy-Green strain tensors 


F = g„ ^ G“, C = [g„ ■ g^j G“ ® G^ C, = [g„ ■ gf,] G“ 0 G 


(33) 


is written in the natural covariant basis as B = 


If a vector 

[ga ■ g^j G“. Thus the identity tensor used for double contraction in (|32|) 


Ga, then CB = 


1,4 


needs to 


have a covariant set of basis vectors while that used in (32)3 requires a contravariant 
basis. 

To obtain the non-equilibrium part of the energy density function, we consider a 


simplihcation of the Mooney-Rivlin type energy in (31) by taking cxe = 0 and n = 1. 


Since the non-equilibrium part of energy should depend only on the elastic parts of the 
deformation gradient and the magnetic induction as assumed earlier, we require R in 
d^i to obtain the value Ce : I. This is equivalent to the expression C : on using 

Eq. ([^. Hence to obtain the non-equilibrium energy, instead of a double contraction with 


the identity tensor to obtain the invariants in (32), we do so by the contravariant tensor 


and the covariant tensor ^ as appropriate and replace the Lagrangian vector B by 
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Bg. The energy function thus obtained is similar to a generalisation of a neo-Hookean 
type energy to include magnetic effects 


n,{C, C„ B, B,) = ^ [c : - 3] + g. 




[C [B - B^] ] 0 [C [B - BJ ] 


: C, 

: C-'. 


(34) 


The magneto-viscoelastic coupling parameters and here are similar to q and r used 
in the dehnition of Vte and we have used the relation Bg = B — B^. 

A strong coupling between the mechanical viscous measure and the magnetic non¬ 
equilibrium quantity B^ is noted here. For the sake of simplicity of our calculations and 
in the absence of any experimental data, we simplify the above expression by assuming 
that the magnetic non-equilibrium effects are coupled only with the total deformation C 
and not with its viscous component C^. This assumption was also used earlier to arrive 


at the dissipation conditions (24). Thus, by replacing with I in (34), a simpler form 


of the non-equilibrium energy is obtained as 


SJ„(C, C„, B, B„) = ^ [C-i : C - 3] + 5, 


.]] 


+rv 


[c 


[c 


J] 


: I, 


(35) 


3.2 Evolution equations 


In order to completely dehne the magneto-viscoelastic behaviour of a solid material, 
along with the balance laws and the energy density functions dehned in the previous 
sub-section, we also require evolution laws for the ‘viscous’ quantities B^ and C^,. These 
are postulated such that the laws of thermodynamics are satished at every instant and 
B.y and stop evolving when the system reaches an equilibrium state. 

For the non-equilibrium part of the magnetic induction we consider the following 


evolution equation such that the left side of inequality (27) becomes a negative semi- 
dehnite quadratic form, automatically satisfying the thermodynamic constraint. Thus, 


dMy Po 

dt Tm 5B„ ’ 

= [94 + »-„C| 

■L m. 


(36) 


For the mechanical viscous strain tensor, we use the evolution equation as proposed 


by Koprowski-Theiss et ah (2011) 


dt 


1 

-L 1) 


C - - [C : cy] C. 


(37) 


In the equations above, Ty is the specihc relaxation time for the viscoelastic component 
of the dissipation mechanism while is the specihc relaxation time for its magnetic 
component. Typically Ty is of the order of some minutes or even upto a few hours while 
Tm is of the order of a few seconds or some milliseconds. 


It now remains to prove that the evolution equation (37) is thermodynamically con 


sistent with the energy density function (35), i.e. they satisfy the constraint (26). 
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Consider a fourth order projection tensor defined as 


Idev 


1 


c 


-1 


where I is the fourth order symmetric identity tensor given in component form as 

\^]ijkl ~ 9 + ^il^jk) i 


(38) 


(39) 


5ij being the Kronecker-Delta. On a double contraction of this tensor with C, a 


multiple of the right side of the evolution law ( |37[ ) is obtained. Thus, 

1 


: C = 


dt 


From Eq. (35), we obtain 

dVt 


sc„ = ■y®’' = y^ ■ 


1-1 


(40) 


(41) 


where we have used the negative-dehnite fourth order projection tensor I^-i which when 
expanded in component form, gives 


[[C.-y, [€„-■],, + [C„-y, [cP],, 


\ijkl 2 

Consider the following operation 

dVLv dC„ jj,- 


ijk 


dC,. 


, = ^ [C : Ic-i] : [Ic" : Cl . 

dt 2T„ L c„ j L c„ j 


Since the tensor is idempotent, the above expression can be rewritten as 

dfl^ dC^ fiyr r .gy 

ac; ■ yr = ^ -F. 

= lip'' : Cl : Ip-i : \lp'' : Cl < 0. 

2^ L J L '“'■u J — 


(42) 

(43) 

(44) 

(45) 


The above inequality holds since Ipj-i is negative dehnite. This is the statement of 


the dissipation condition (26), hence the evolution equation (37) is thermodynamically 


consistent with the energy density function (35). 


3.3 Stress and magnetic field calculations 


For the energy functions dehned in Eqs. (31) and (35), the Piola-Kirchhoff stress is given 
as 


S = 2 




+ 2 


,dilv 


where 


Se = ^ 
2 


1 + cxe tanh ( — 


pC-^ = Se + S„-pC-^ 

[[l + n]I+[l-n] [/iI-C]] 
+2rB ® [CB] + 2r [CB] (g) B, 


(46) 


(47) 
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and 


Si; — fiyC^ + 2ri,]Be ® [CBe] “1“ 2ri, [CBe] 
The Lagrangian magnetic field El is given as 


e = 


dflf. ^ dflr 


aB dM 


= Elg + El,;, 


(48) 


(49) 


where 


H, = 


he 

2me 


I — tte tanh I — 
rrie 


[[1 + n] [/,-3] + [l-n] [12-3]]] 
+2# + 2rC^I 


and 


(50) 


(51) 


Ell, = 2g„Be + 2r„CBe. 

The ‘viscons’ or non-eqnilibrinm magnetic field defined above tends to zero at eqnilibrinm 
when Be —)■ 0. Enlerian expressions for the eqnilibrinm valnes of the total Canchy stress 


T and the magnetic field E can be written nsing Eqs. (47) and (50) as 


he 

Te = — 


1 + tte tanh I — 

rrip 


[1 + n] b + [1 — n] [Jib — b^] 
+r b ® [bb] + r [bb] ® b, 


(52) 


bp = 


cxetanh"' ( — 

' , mp 


[ [1 + n] [Ji - 3] 


he 
2me 

+ [1 — n] [I 2 — 3] ] b“^b + 2gb“^b + 2rbb. 


(53) 


In the case of no deformation (b = i), if r = 0, then the total eqnilibrinm stress in 


Eq. (52) is nnaffected by the magnetic indnction and if g = r = 0 then the eqnilibrinm 


magnetic field in Eq. (53) is nnaffected by the nnderlying deformation. The conpling 


cansed by the parameters q and r between deformation and magnetic field is inverse to 
each other. The former canses a directly proportional relation of b~^ to Eg while the 
latter links b to Eg. Similarly and r„ are reqnired for inclnding these two-way conpling 
effects for the non-eqnilibrinm qnantities. 

For the case of no deformation of an isotropic material, the magnetic field shonld be 
in the direction of the applied magnetic indnction and be directly proportional to the 


latter. From Eq. (53), this imposes the constraint 

g -t- r > 0. 


(54) 


If r > 0, the material stiffens in the direction of the applied magnetic indnction while 
if Og > 0, the total stiffness of the material increases isotropically. Both these effects have 


been observed to be the case in many MREs, see, for example, the resnlts of Jolly et ah 


(1996 

) and 

Varga et al.j ( 

2006 


However, this need not necessarily be trne in general 
for all magnetoelastic materials and r can have negative valnes. In this sitnation and 
the case of a material with weak magnetoelastic conpling, i.e. very small valnes of r; we 


reqnire g > 0 to satisfy the constraint (54). 


12 


























4 Numerical examples 

In this section, we model four different types of experiments and obtain the corresponding 
solutions numerically. The following numerical values of the material parameters are used 
unless otherwise stated to have a different value for individual computations 

yUo = dvr X 10“'^ N/A^, /Xg = 2.6 x 10^ N/m^, /x^ = 5 x 10® N/m^, 

ag = 0.3, mg = lT^ n = 0.3, g = r = = l//xo, = 5//io. (55) 


The value of /Xg is taken to be the value of shear modulus at zero magnetic held for an 


elastomer hlled with 10% by volume of iron particles, cf. 

Jolly et ah 

(1996 

). Values 

of n, g and r are what have been used by 

Ottenio et ah 

(2008 

) and 

Saxena and Ogden 


(2011). Values of /x„, Ofg, g^, are within reasonable physical assumptions and we analyse 
the dependence of our solutions on the values of these parameters. 

For computations in the following subsections, the equations derived earlier are spe¬ 
cialised to a uniaxial deformation and magnetisation in cartesian coordinates. The time- 
integration is performed using a standard solver ode45 from Matlab that employs an 
explicit Runge-Kutta scheme, cf. Shampine and Reichelt (1997). 


4.1 Magnetic induction with no deformation 


With a motivation to isolate and understand the effects of the applied magnetic induction 
on the magneto-viscoelastic deformation process, we consider no deformation in this hrst 
case. Consider an experiment with the sample held hxed at zero deformation (Ai = 
-^2 = -^3 = 1) and a sudden but constant magnetic induction applied at time t = 0. 
This results in the generation of a viscous overstress and a temporary increment in the 
magnetic held, both of which settle down to equilibrium values with time. Variations of 
the total magnetic held Ihi (component of Ih in the xi direction) and the total Cauchy 
stress Til with time are plotted in Figs. [2]-[^ for the following values of the magnetic 
induction. 


®2 — ®3 — 0 , 




0, for f < 0, 
0.1 T, forf>0. 


(56) 


We study the dependence of the magneto-viscoelastic coupling parameters g^ and 
r^,, and the applied magnetic induction Bi on the relaxation of magnetic held and total 
Cauchy stress. 

It is seen from Fig. that a large value of causes a high initial magnetic held 
but the decay to equilibrium value is also faster when g„ is high. In the case of the total 
Cauchy stress, as seen from Fig. [^, g^ has no ehect on the initial viscous overstress but 
a large g^ also causes the stress to decay faster and reach the equilibrium value. This is 
expected since the expression for in Eq. (48) does not contain g^ explicitly but the 
dependence comes through the evolution equation (36). The parameter has a similar 
ehect on the magnetic held as does but diherent in the case of the total Cauchy stress. 
As observed from Fig. a large value of causes a higher initial stress and a faster 
decay of the same to equilibrium. The dependence of the relaxation processes on the 
applied magnetic induction is shown in Fig. It is seen that the equilibrium values for 
all curves are diherent in this case since they depend on the value of applied magnetic 
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Fig. 2. Variation of (a) the total magnetic field Ihi (A/m) and (b) the principal total 
Cauchy stress rn (N/m^) with time t (s) for no deformation in the presence of a step 
magnetic induction Bi = 0.1 T. Four curves correspond to different values of (i) = 

1/ho, (ii) qv = 2//io, (hi) = 4//io, (iv) q^ = 7/po- 




Fig. 3. Variation of (a) the total magnetic held Ihi (A/m) and (b) the principal total 
Cauchy stress tu (N/m^) with time t (s) for no deformation in the presence of a step 
magnetic induction Bi = 0.1 T. Four curves correspond to different values of (i) = 

1/po, (h) Ty = 2//io, (hi) ry = 4/po, (iv) Vy = T/po- 
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Til 



Fig. 4. Variation of (a) the total magnetic field Ihi (A/m) and (b) the principal total 
Cauchy stress th (N/m^) with time t (s) for different values of the magnetic induction 
Bi (i) Bi = 0.05 T, (ii) Bi = 0.1 T, (hi) Bi = 0.2 T, (iv) Bi = 0.3 T. 


induction. As expected from Eqs. (52) and (53) a higher magnetic induction causes a 
larger magnetic held and a larger stress. It is also observed that the higher the magnetic 
induction, the longer it takes for both the magnetic held and the stress to relax and reach 
equilibrium. 


4.2 Magnetic induction with a uniaxial deformation 

In this case, we specify a deformation and a magnetic induction in the xi direction while 
allowing the material to move freely in X 2 and x^ directions. The stretch Ai and the 
magnetic induction Bi are applied at time t = 0 and then the material is allowed to relax 
and reach an equilibrium state. Variation of the magnetic held and the total Cauchy 
stress with time are plotted in Figs. and for the values 

r^ = 100s, Bi = 0.1T, Ai = 1.5, (57) 


unless otherwise stated to be diherent for individual problems. 

It is observed from Fig. that two diherent relaxations over diherent time scales 
occur in the total Cauchy stress Th ~ one corresponding to the evolution of B^ and other 
corresponding to that of C^,. The decay for small time scale (upto 1.5 seconds) to reach 
an equilibrium value is shown in Fig. while that for a longer time scale (upto 1000 
seconds) is shown in Fig. |^. The four curves correspond to four diherent values of the 
initial stretch. It should be noted that the end point of a curve in Fig. is the same 
as the starting point of the corresponding curve in Fig. [^. A higher stretch causes an 
increase in the equilibrium value of the magnetic held in Fig. and also causes a faster 
relaxation to equilibrium value. As expected from the existing results of pure mechanical 
viscoelasticity, cf. Hossain et al. (2012), a higher stretch leads to a larger value of stress 
in Fig. [^and a smaller value of causes stress to relax faster in Fig. |^. 
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Fig. 5. Uniaxial stretch in xi direction, Bi = 0.1 T. Variation of principal total Cauchy 
stress Til (N/m^) vs time t (s) for different values of the stretch Ai. (a) small time scale 
(b) large time scale, (i) Ai = 1.5, (ii) Ai = 2, (hi) Ai = 3, (iv) Ai = 4. 




Fig. 6. Uniaxial stretch in xi direction, Bi = 0.1 T. Variation of (a) total magnetic held 
111 (A/m) vs time t (s) for different values of the stretch Ai. (i) Ai = 1.5, (ii) Ai = 2, (iii) 
Ai = 3, (iv) Ai = 4; (b) principal total Cauchy stress th (N/m^) vs time t (s) for different 
values of parameter (i) = 50 s, (ii) = 100 s, (iii) = 200 s, (iv) T„ = 300 s. 
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Fig. 7. Principal total Cauchy stress rn (N/m^) vs stretch Ai. (a) Different values of 
the underlying magnetic induction (i) Bi = 0, (ii) Bi = 0.1 T, (iii) Bi = 0.3 T, (iv) Bi = 
0.5 T. (b) Different values of the stretch rate (i) Ai = ±0.005 s“^, (ii) Ai = ±0.02 s“^, 
(iii) Ai = ±0.04 s“^, (iv) Ai = ±0.08 s“^. 


4.3 Time dependent deformation 

We now study the effects of the magneto-viscoelastic coupling on a dynamic deformation 
of the material. In this case, calculations are performed corresponding to an experiment 
where a magnetic induction is applied at time t = 0 and the material is stretched with 
a constant rate in the xi direction. On reaching Ai = 3, stretch is reduced at the same 
rate until a condition of zero stress or zero deformation (whichever earlier) is reached. 
Effects on the total Cauchy stress and the total magnetic field of the applied magnetic 
induction, the rate of stretch, and the parameters and is analysed in Figs. 01 The 
following values of the magnetic induction and the stretch rate are used 

Bi = 0.2T, Ai = ±0.01s“\ (58) 


unless otherwise stated to be different for individual calculations. 

It is seen from Fig. that the starting points of all four curves are different corre¬ 
sponding to the stress induced due to the applied magnetic induction. The stress first 
increases with time (due to an increasing Ai) and then falls with a decreasing Ai following 
a different path than earlier. A higher magnetic induction leads to larger value of the 
peak stress reached during the process. Similar curves for different values of the stretch 
rates are shown in Fig. [^. A larger value of stretch rate causes a larger peak value of 
stress since the material gets less time to relax as observed for the purely mechanical 
viscoelastic case by Lion (1997) and Amin et ah (2006). 

Similar variation of the magnetic field Ihi with the stretch Ai can be observed for large 
values of stretch rates since is much smaller than T^. The results for these calculations 
are shown in Figs. and for Ai = 2 as the value of the maximum obtained stretch. As 
observed from Fig. starting at t = 0, the magnetic field first falls and then rises due 
to an increase in Ai. As Ai reduces, Ihi comes down approaching a steady equilibrium 
value. High values of and cause a high initial magnetic field and a faster approach 
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Fig. 8. Total magnetic field Ihi (A/m) vs stretch Ai at a stretch rate Ai = ±3 s~^. 
(a) Different values of the parameter (i) = l/fio, (ii) q^ = 3//io, (hi) qv = 5//io, 

(iv) q^ = lO/yUo; (b) Different values of the parameter (i) = l/po, (u) = 3//io, 

(hi) Ty = 5//io, (iv) Ty = 10//io. 




Fig. 9. Total magnetic field Ihi (A/m) vs stretch Ai. (a) Different values of the magnetic 
induction Bi at a stretch rate Ai = ±3 s“^ (i) Bi = 0.1 T, (ii) Bi = 0.2 T, (hi) Bi = 0.3 
T, (iv) Bi = 0.4 T; (b) Different values of the stretch rate Ai (i) Ai = 1 s“\ (ii) Ai = 2 
s“^, (hi) Ai = 3 s“^, (iv) Ai = 4 s“^. 
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Fig. 10. Time dependent magnetic induction, variation with the rate of magnetic induc¬ 
tion (a) Total magnetic held Ihi (A/m) vs magnetic induction Bi (T) (i) Bi = 1 T/s, 
(ii) Bi = 2 T/s, (iii) Bi = 3 T/s, (iv) Bi = 4 T/s; (b) Principal total cauchy stress 
Til (N/m^) vs magnetic induction Bi (T). (i) Bi = 1 T/s, (ii) Bi = 4 T/s. 


towards the equilibrium. 

Dependence of this process on the applied magnetic induction and the stretch rate is 
shown in Fig. The different start and end points of the curves in Fig. correspond to 
the values of magnetic held caused by diherent magnetic inductions. A higher magnetic 
induction causes a larger magnetic held while for a lower value of stretch rate, as observed 
from Fig. |^, the magnetic held approaches the steady equilibrium value earlier in the 
cycle. 


4.4 Time dependent magnetic induction 


In this case we study the ehect of a time-varying magnetic induction on the induced stress 
and magnetic held in an undeformed material. The sample is assumed to be hxed at zero 
deformation and a magnetic induction is applied at time f = 0 in the Xi direction with 
a constant rate until a value of Bi = 0.8 T is obtained. The induction is then reduced 
with the same rate until it reaches zero. Numerical results for this case are shown in 
Figs. 10 -13 A value of Bi = 2 T/s is used to plot the curves in Figs, ll]- 13 


It is observed from Fig. [TO^ that for a particular value of induction rate, the magnetic 
held increases with an increasing magnetic induction and in the return cycle, it reduces 
to eventually obtain a negative value. In the entire cycle, the material develops a mag¬ 
netisation in the xi direction due to the existing magnetic induction. As the induction 
reduces to zero, the material has to develop a magnetic held in the negative xi direction 
to erase the magnetisation in xi direction. Also, it can be seen by Eq. (|^ that if b = 0 , 
then Ih and m obtain opposite signs. 

A higher rate of magnetic induction causes the magnetic held to reach a high peak 
value and it also causes the magnetic held to reach the maximum negative value when 
B vanishes completely. Stress in this case has a rather interesting variation as observed 
from Fig. [T0|3. Starting from zero, the stress increases with an increasing Bi and then 
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Fig. 11. Time dependent magnetic induction, variation with the parameter (a) Total 
magnetic held Ihi (A/m) vs magnetic induction Bi (T) (i) = 0.1//io, (h) qv = l//^o, 

(iii) = 5//io; (b) Principal total cauchy stress rn (N/m^) vs magnetic induction Bi (T). 
(i) q^ = 0.1//io, (ii) qv = ^/no- 




(b) 


IBi 


Fig. 12. Time dependent magnetic induction, variation with the specihc relaxation time 
Tm (a) Total magnetic held Ihi (A/m) vs magnetic induction Bi (T) (i) = 1 s, 

(ii) Tm = 2 s, (iii) Tm = 3 s; (b) Principal total cauchy stress rn (N/m^) vs magnetic 
induction Bi (T). (i) Tm = 1 s, (ii) Tm = 3 s. 
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Fig. 13. Variation of the total magnetic field Ihi (A/m) with the magnetic induction 
Bi (T), Bi = 2 T/s, (i) Ai = 1, (ii) Ai = 1.2, (iii) Ai = 1.4, (iv) Ai = 1.6. 


falls as Bi is reduced to zero. However, for a high induction rate of Bi = 4 T/s, in the 
return cycle the stress approaches a minimum and then rises again. 

A larger value of causes a higher value of peak magnetic held and a larger negative 
value at the end of cycle as Bi —)■ 0. Moreover, it also clearly causes larger energy 
dissipation during the cycle as the area inside the curve (iii) of Fig. Hr is much larger 
than that in curve (i). The stress in Fig. 11 d has a higher peak value for a smaller value 
while a larger value of 


of g. 

value of the parameter T, 


causes a higher stress at the end of the cycle. A smaller 
helps the material to relax in lesser amount of time, hence the 


peak values of the magnetic held and stress reached in Fig. 12 in this case are lower. A 
higher stretch in Fig. 13 causes a larger peak value of the magnetic held. 


5 Concluding remarks 

We have presented a theory to model nonlinear magneto-viscoelastic deformations in 
this paper. The deformation gradient is multiplicatively decomposed and the magnetic 
induction is additively decomposed to ‘elastic’ and ‘viscous’ parts to take into account 
dissipation mechanisms. A Mooney-Rivlin type magnetoelastic energy density function 
is used for the equilibrium part, which is simplihed to a neo-Hookean type energy density 
function for the non-equilibrium part of the free energy. These, along with thermody¬ 
namically consistent evolution laws, are used to obtain numerical solutions corresonding 
to several diherent magneto-viscoelastic deformations. 

The magneto-viscoelastic parameters q^ and can have strong ehects on the non- 
equilibrium magnetic held and the non-equilibrium total Cauchy stress by changing their 
peak values and the decay times. Strong couplings are also shown to exist between the 
magnetic induction and the non-equilibrium stress, and the underlying deformation and 
the non-equilibrium magnetic held, as is evident from Figs. 4b and 6a. We observe that 
a stretch rate and a magnetic induction rate can have a considerable inhuence on the 
total Cauchy stress and the magnetic held. The developed model seems to capture the 
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magneto-viscoelastic phenomena quite nicely and on isolating mechanical viscoelastic 


effects, our results are qualitatively the same as those obtained earlier by Amin et ah 
(2006) and Hossain et ah (2012). 

It should be noted that the numerical results presented here are representative so¬ 
lutions considering only one dissipation mechanism in the body. The theory can be 
easily generalised to include multiple mechanisms to match the experimental data. We 
have considered a specific type of phenomenologically motivated constitutive law for an 
isotropic material in this paper. The theory for an anistropic material and possibility of 
existence of other constitutive laws (such as those derived from micromechanics of the 
material) will be discussed in forthcoming contributions. 
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